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ABSTRACT 

We study the leading corrections to the emergent canonical com- 
mutation relations arising in the statistical mechanics of matrix models, by 
deriving several related Ward identities, and give conditions for these cor- 
rections to be small. We show that emergent canonical commutators are 
possible only in matrix models in complex Hilbert space for which the num- 
bers of fermionic and bosonic fundamental degrees of freedom are equal, 
suggesting that supersymmetry will play a crucial role. Our results simplify, 
and sharpen, those obtained earlier by Adler and Millard. 
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1. Introduction 

It is widely believed that at distances of order the Planck length £p ~ 10~^^ cm our 
conventional notions of the geometry of spacetime break down, as a result of quantum gravity 
effects. One indication of the modifications in physics that might be expected is provided 
by string theory models of quantum gravity, in which several studies suggest a modification 
of the uncertainty relation of the form [1] 

AxAp>^[l + P{Apf + ....], p>0 , (la) 

implying a finite minimum uncertainty Axq — hp^ in the vicinity of the Planck length. 
As discussed by Kempf and collaborators [2], Eq. (la) corresponds to a correction to the 
Heisenberg canonical commutation relations of the form 

[x,p] = in{i + (3p^ + ...) . (lb) 

We wish in this paper to discuss modifications of the Heisenberg algebra arising in 
another context, that of the statistical mechanics of matrix models, and to compare them 
with Eqs. (la, b). Several years ago, Adler proposed a set of rules for a generalized quantum 
or trace dynamics, which is a Lagrangian and Hamiltonian mechanics with arbitrary non- 
commutative phase space variables q,p, and this was developed in a series of papers with 
various collaborators [3]. For theories in which the action is constructed as the trace of a 
sum of matrix products of N x N matrix variables, trace dynamics gives a powerful, basis 
independent, way of representing the same dynamics that can also be described in terms 
of the N'^ individual matrix elements. A significant new result emerging from this point of 
view was obtained by Adler and Millard [4], who argued that the statistical mechanics of 
trace dynamics takes the form of conventional quantum field theory, with the Heisenberg 
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commutation relations holding for statistical averages over certain effective canonical vari- 
ables obtained by projection from the original operator canonical variables. Recently, it has 
become clear [5] that the underlying assumptions of trace dynamics are satisfied by matrix 
models, for which the methods of trace dynamics provide a very convenient calculational 
tool. Hence the results of Adler and Millard can be reinterpreted as providing a statistical 
mechanics of matrix models, and showing that thermal averages in this statistical mechanics 
can behave as Wightman functions in an emergent local quantum field theory. These results, 
together with recent work [6] suggesting that the underlying dynamics for string theory may 
be a form of matrix model, raise in turn the question of determining the form of the leading 
corrections to the Heisenberg algebra implied by the statistics of matrix models, formulating 
conditions for these corrections to be small, and seeing whether they can be related to the 
string theory result of Eqs. (la, lb). 

An investigation of these questions is the focus of this paper, which is organized as 
follows. In Sec. 2 we give a brief synopsis of the rules of trace dynamics in the context of 
matrix models. We show that the conservation of the operator [4, 7] C can be understood as 
a simple consequence of unitary invariance. We also remark that, with Grassmann fermions, 
C is independent of the classical parts of the matrix phase space variables, and review the 
statistical mechanics [4, 8] of matrix models. In Sec. 3 we consider the simple case of a 
bosonic matrix model with Hamiltonian quadratic in the canonical momenta, and, making 
no approximations, derive a simplified form of the Ward identity used in Ref. [4] to obtain the 
effective canonical algebra. This analysis shows that there are corrections to the canonical 
commutator quadratic in the canonical momentum. In Sec. 4 we use the symplectic formalism 
of Ref. [4] to repeat this calculation in the case of a general commutator/anticommutator of 
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canonical variables in a generic matrix model, that can include fermions, and we generalize 
the treatment of [4] to allow nonzero sources for the classical parts of the matrix variables. 
Prom the analyses of Sees. 3 and 4, we formulate conditions for the corrections to the 
emergent canonical algebra to be small. We show that these conditions require C to be 
an intensive rather than extensive thermodynamic quantity, and that they can be satisfied 
in complex Hilbert space (if at all) only in matrix models with precisely equal numbers of 
bosonic and fermionic degrees of freedom. This result strongly suggests that candidate matrix 
models for prequantum mechanics should be supersymmetric. We conclude by generalizing 
the conditions to ones that permit the recovery of the full emergent quantum field theory 
structure derived in Ref . [4] . We also compare the prequantum corrections to the canonical 
algebra derived in Sees. 3 and 4, in which (as is usual in field theories) the spatial coordinate 
is simply a label, and the field variables are the dynamical canonical variables, to the string 
theory inspired expression of Eq. (lb) in which a; is a coordinate operator. 

2. The Statistical Mechanics of Matrix Models 

We begin by reviewing the statistical mechanics of trace dynamics, taking into ac- 
count the simplifications [5] that become possible when Grassmann algebras are employed 
to represent the fermion/boson distinction. Let Bi and B2 be two N x N matrices with 
matrix elements that are even grade elements of a complex Grassmann algebra, and Tr the 
ordinary matrix trace, which obeys the cyclic property 

TrBiB2 = ^(-Bl)mn(-B2)nm = ^(-B2)nm(-Bl)mn = TrS2-Bi . (2a) 

m,n m,n 

Similarly, let xi X2 be two N x N matrices with matrix elements that are odd grade 
elements of a complex Grassmann algebra, which anticommute rather than commute, so that 
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the cyclic property for these takes the form 

TrXlX2 = 5Z(Xl)mn(X2)nm = " ^{X2) nm(Xl) mn = -T^X2Xl ■ (2&) 
m,n m,n 

The cychc/ anticychc properties of Eqs. (2a, 2b) are just those assumed for the trace operation 
Tr of trace dynamics.* From Eqs. (2a, b), one immediately derives the trilinear cyclic 
identities 

TrBi[B2, S3] =IVS2[S3, B,] = TrBs[Bi, B2] 
TrSi{S2, 53} =TrB2{B^, B^} = TrB^iB^, B^] 

(2c) 

Tr5{xi,X2} =Trxi[x2,5] = Trx2[Xi,^] 

IVxi{S,X2}=TV{xi,5}X2 = TV[xi,X2]S , 
which are used repeatedly in trace dynamics calculations. 

The basic observation of trace dynamics is that given the trace of a polynomial P 
constructed from noncommuting matrix or operator variables (we shall use the terms "ma- 
trix" and "operator" interchangeably in the following discussion) , one can define a derivative 
of the c-number TrP with respect to an operator variable O by varying and then cyclically 
permuting so that in each term the factor 50 stands on the right, giving the fundamental 
definition 

(5TrP = TX——50 , (3a) 
60 



* In Refs. [3, 4] the fermionic operators were realized as ordinary matrices with com- 
plex matrix elements, all of which anticommute with a grading operator (—1)^ which formed 
part of the definition of the graded trace Tr, for which fermions then obeyed Eq. (2b) while 
bosons obeyed Eq. (2a). Since the use of Grassmann odd fermions eliminates the need for 
the inclusion of the (—1)^ factor, the graded trace obeying Eqs. (2a, b) is here just the usual 
matrix trace, for which we use the customary notation Tr. 



or in the condensed notation that we shall use throughout this paper, in which P = TrP, 

SP^Tr^^50 . {3b) 

Letting L[{gr}, {<ir}] be a trace Lagrangian that is a function of the bosonic or fermionic 
operators {g^} and their time derivatives (which are all assumed to obey the cyclic relations 
of Eqs. (2a-c) under the trace), and requiring that the trace action S = J dtL be stationary 
with respect to variations of the q^s that preserve their bosonic or fermionic type, one finds 
[3] the operator Euler-Lagrange equations 

Sh d 6L ^ 

5qr dt 6qr 

Because, by the definition of Eq. (3a), we have 

r6L\ _ dL 
\SqrJij~ d{qr)ji 

for each r the single Euler-Lagrange equation of Eq. (3c) is equivalent to the N'^ Euler- 
Lagrange equations obtained by regarding L as a function of the iV^ matrix element variables 
{qr)ji- Defining the momentum operator pr conjugate to qr, which is of the same bosonic or 
fermionic type as qr, by 

Pr = T^ : (4a) 

dqr 

the trace Hamiltonian H is defined by 

U^TrJ^PrQr-'L . (4&) 

r 

In correspondence with Eq. (3d), the matrix elements {pr)ij of the momentum operator pr 
just correspond to the momenta canonical to the matrix element variables {qr)ji- Performing 



general same-type operator variations, and using the Euler-Lagrange equations, we find from 
Eq. (4b) that the trace Hamiltonian H is a trace functional of the operators {qr} and {pr}, 

H = H[{g,},{p,}] , (5a) 

with the operator derivatives 

with er = 1(— 1) according to whether qr,Pr sue bosonic (fermionic). Letting A and B be 
two trace functions of the operators {g^} and {pr}, it is convenient to define the generalized 
Poisson bracket 

r. X.. .T. V- fSASB 5B5A\ 

{A,B} = TrV eJ — ■ 6a 

\dqr OPr dqrOPrJ 

Then using the Hamiltonian form of the equations of motion, one readily finds that for a 
general trace functional A[{gr}, {pr}], the time derivative is given by 

|a = {A,H} ; (66) 

in particular, letting A be the trace Hamiltonian H, and using the fact that the generalized 
Poisson bracket is antisymmetric in its arguments, it follows that the time derivative of H 
vanishes. An important property of the generalized Poisson bracket is that it satisfies [3] the 
Jacobi identity, 

{A,{B,C}} + {C,{A,B}} + {B,{C,A}} = . (6c) 
As a consequence, if Qi and Q2 are two conserved charges, that is if 

0=^Qi = {Qi,H} , 0=^Q2 = {Q2,H}, (6d) 



then their generahzed Poisson bracket {Qi,Q2} also has a vanishing generahzed Poisson 
bracket with H, and is conserved. This has the consequence that Lie algebras of symmetries 
can be represented as Lie algebras of trace functions under the generalized Poisson bracket 
operation. 

A significant feature of trace dynamics is that, as discovered by Millard [7], the 
anti-self-adjoint operator [7, 4] 

C= ^ [qr,Pr\- XI {<lr,Pr} (7) 
r bosons r fermions 

is conserved by the dynamics. As we shall now show, conservation of C holds whenever 
the trace dynamics has a global unitary invariance, that is, whenever the trace Hamiltonian 
obeys 

U.[{UUrU}, {V^PrU}] = H[{?,}, {pr}] (8a) 

for a constant unitary N x N matrix [/, or equivalently, by Eq. (4b), whenever the trace 
Lagrangian obeys 

L[{C/tg,C/}, {UHrU}] = L[{g,}, {g,}] . (86) 

Letting U — exp A, with A an anti-self- adjoint bosonic generator matrix, and expanding to 
first order in A, Eq. (8a) implies that 

H[{g, - [A, Qr]}, {pr - [KPr]}] = H[{g,}, {pr}] . (9a) 

But applying the definition of the variation of a trace functional given in Eq. (3b), Eq. (9a) 
becomes 

T'E(-fiA.*i-f|A.P'i)=° . m 
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which by use of the trihnear cychc identities of Eq. (2c) yields 

TrA> -^gr-Cr^rT h "^Pr " CrPr"^ = 0. (9c) 

^ \dqr OQr OPr OPr J 

Since the generator A is an arbitrary anti-self-adjoint N x N matrix, the anti-self- adjoint 
matrix that multiplies it in Eq. (9c) must vanish, giving the matrix identity 

/ ^ '^Ir - ^rQr-r V -^—Pr " (^rPvT— = 0. (lOo) 

\Oqr OQr OPr OPr J 

But now substituting the Hamilton equations of Eq. (5b), Eq. (10a) takes the form 

= J]] {-Prqr + erQrPr + ^rQrPr - PAr) 
r 

— -^ 5^ (—PrQr + ^rQrPr) (106) 



( ilr^Pr]- {^r,Pr}] , 

\r bosons r fermions / 



completing the demonstration of the conservation of C. 

Corresponding to the fact that C is conserved in any matrix model with a global 
unitary invariance, it is easy to see [4, 8] that C can be used to construct the generator of 
global unitary transformations of the Hilbert space basis. Consider the trace functional 

Ga = -TrAC , (11a) 

with A a fixed bosonic anti-self-adjoint operator, which can be rewritten, using cyclic invari- 
ance of the trace, as 

GA = Tr^[A,p,]g, = -Tr2];p,[A,g,] . (116) 

r r 

Hence for the variations of Pr and Qr induced by Ga as canonical generator, which have a 
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structure analogous to the Hamilton equations of Eq. (5b), we get 



sol 

6qr=er^ — = -[A,qr] . 



Comparing with Eqs. (8a) and (9a), we see that these have just the form of an infinitesimal 

global unitary transformation. 

For each phase space variable qr,Pr, let us define the classical part q'^-iP^ and the 

noncommutative remainder q'^ip'^, by 

1 1 

^ ^ (12a) 

Qr ^Qr - qr Pr ^ Pr - Pt , 

so that bosonic q^jP^ are c-numbers, fermionic q^jP^ are Grassmann c-numbers (where by 
a c-number we mean a multiple of Ijv, the N x N unit matrix), and the remainders are 
traceless, 

Trg; = Trp; = . (126) 

Then since q^jPr commute (anticommute) with qg,p'g for r, s both bosonic (fermionic), we see 
that the classical parts of the phase space variables make no contribution to (7, and Eq. (7) 
can be rewritten as 

C= Yl ^r^P'r]- E i^r^P'r} " (12^) 

r bosons r femions 

Thus C is completely independent of the values of the classical parts of the matrix phase 
space variables. 

Making the assumption that trace dynamics is ergodic (which undoubtedly requires 
an interacting as opposed to a free theory, and may presuppose taking the iV — > oo limit), one 
can then analyze [4] the statistical mechanics of trace dynamics for the generic case in which 
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the conserved quantities are the trace Hamiltonian H and the operator C. As discussed in 
detail in the second paper cited in Ref. [5], the analysis of [4] carries over to the case in which 
the fermions are represented by Grassmann matrices; the demonstration of a generalized 
Liouville theorem still holds, and the requirements for convergence of the partition function 
are much less stringent, eliminating the complexities addressed in Appendix F of Ref. [4]. 
With Grassmann fermions, for the typical models we are studying the bosonic part of i7 is a 
positive operator, from which H inherits good positivity properties . The canonical ensemble 
then takes the simple form given in Eq. (48c) of [4] , 



with dii the invariant matrix (or operator) phase space measure provided by Liouville's the- 
orem, with T a real number, and with A an anti-self-adjoint matrix that in the generic case 
(which we assume) has no zero eigenvalues. (Equation (13) can be derived directly [4] by 



or indirectly [8] by first calculating the corresponding microcanonical ensemble correspond- 
ing to these conserved quantities, and then using standard statistical physics methods to 
calculate the canonical ensemble from the microcanonical one.) We wish to make two points 
about the partition function defined in Eq. (13). First of all, it is not invariant under the 
unitary transformation of Eq. (8a) for fixed A, but is invariant when A is simultaneously 
transformed to WXU; hence the partition function breaks unitary invariance, but has a 
specific form of unitary covariance. Second, the partition function contains a weighted sum 
over all possible commutators [q'r,Ps] for bosonic variables and all possible anticommutators 
{qr,Ps} for fermionic variables; there is no restriction to the classical or quantum mechanical 
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p=Z-^exp{-TH-TrXC) 




(13) 



maximizing the entropy subject to the constraints imposed by the conservation of H and C, 



evaluation of these commutators/anticommutators as or i5rs respectively. However, sta- 
tistical integrals like Eq. (13) are typically dominated by specific regions of the integration 
domain, and we will see, by a study of the Ward identities following from Eq. (13), that 
this can lead to effective quantum mechanical commutators inside statistical averages. The 
structure of the Ward identities or equipartition theorems following from Eq. (13) will be 
reanalyzed in the next two sections without making approximations used in [4], so as to 
determine the leading corrections to the emergent canonical commutation relations. Prom 
this analysis we will infer a set of conditions for obtaining the full emergent quantum field 
theory structure of [4]. 



13 



3. Corrections to the Bosonic Commutator [gs,Pr] in 
a Simplified Unitary Invariant Matrix Model 

We consider in this section the simplified bosonic matrix model with trace Hamilto- 



nian 



H = Tr 



(14a) 



2' 

r 

with the qr self-adjoint NxN complex matrix variables and with V a global unitary invariant 
potential. This form is general enough to include the matrix model forms of the bosonic field 
theories of greatest interest, including the Goldstone model, non-Abelian gauge models, and 
the Higgs model. As we saw in the previous section, the Hamiltonian dynamics for this model 
conserves both the real number H and the matrix C, which in this case is given simply by 

C^J2i^r,Pr] . (146) 

r 

Letting p and Z be respectively the canonical ensemble and partition function given in terms 
of H and C by Eq. (13), we define the ensemble average of an arbitrary function O of the 
dynamical variables by 

{0)av = J d/ipO = J d/ie-^^-^^'^^^O . (15) 

Letting O be the conserved operator C, and noting that the right hand side of Eq. (13) can 
be a function only of the ensemble parameters A and r, we have 

{C)Av = fCKr), (16a) 

with / an anti-sclf-adjoint matrix, which in general can be written as a phase matrix i^g 
times a commuting magnitude matrix |/|, 

f^ies\f\, iefr = -l, is^-ieS, [ieff, |/| ] = . (16&) 
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We shall now specialize to an ensemble for which the magnitude matrix |/| makes 
no distinction among the different bases in Hilbert space, and so takes the form of a positive 
real multiple (which we shall call h) of the unit matrix. [As discussed in Appendix B of Ref. 
[4], when H can be expressed in terms of the phase space operators {qr,Pr} using only real 
number coefficients, this assumption implies that we are restricting attention to the special 
class of ensembles for which A = iesA, with A a real multiple of the unit matrix.] Equations 
(16a, b) then become 

{C)av = iesh ■ (16c) 

Since for finite N we necessarily have TrC = 0, the phase matrix must have vanishing 
trace, 

IVieff^O , (16d) 

which implies that ies has N/2 eigenvalues i and N/2 eigenvalues —i. Thus, we are making 
a choice of ensemble for which the U{N) symmetry of H is broken, by the term TrAC, 
to U{N/2) X U{N/2) X R, with R the discrete reflection symmetry that interchanges the 
eigenvalues ±i of iefr- This is clearly the largest symmetry group of the ensemble for which 
one can have {C)av Oj if one were to attempt to preserve the full U{N) symmetry by 
taking an ensemble with A = iA, with A a c-number, then in the canonical ensemble the term 
TrAC would vanish by virtue of the tracelessness of C, and the resulting ensemble would 
have {C)av — 0. Requiring the largest possible nontrivial symmetry group plays the role 
in our derivation of giving a single Planck constant for all pairs of canonical variables; if on 
the other hand, we were to sacriflce all of the U (N) symmetry by allowing generic A, then 
the emergent canonical commutation relations derived below would generically yield N/2 
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different ^'s for the N/2 pairs of canonical variables. It would clearly be desirable to have a 
deeper justification from first principles of our choice of ensemble, perhaps based on a more 
detailed understanding of the underlying dynamics, but at present we must simply introduce 
it as a postulate. 

For this choice of ensemble, let us now consider the Ward identity obtained from 

Z{TYCpr)Av = J dfie-^^-'^~^^TrCpr, (17a) 

by using invariance of the measure d/i under a constant shift of Ps, which implies 
Q=J diidp^ [e-^^-^^'^TrC'p, 

d/xe--"-^^^ (176) 

X [(-r5p,H - Tr~X5p^C)TrCpr + Tr{5p^C)pr + TrC5rsSps] ■ 
Now from Eqs. (14a, b) we have 

Sp^U = Trpsdps, Sp^C = [qs,dps] . (18) 

Substituting these into Eq. (17b), multiplying by Z~^, and using the trilinear cyclic identities 
of Eq. (2c), we get 

= {{-rTrpJp, - Tr[X, q,]5ps)TrCpr + Tr[pr, QsjSps + TrC5rsSps)AV , (19a) 

which since Sps is an arbitrary self-adjoint matrix, implies that the operator multiplying 5ps 
inside the trace must vanish, 

= {{-TPs - [A, qs])TvCpr + \pr: Qs] + C5rs)AV ■ (19^) 

Since A is a constant matrix, it can be taken outside the ensemble average, and so the second 
term in Eq. (19b) takes the form 

-[\{qsTTCpr)Av] , (19c) 
16 



which vanishes since the ensemble average inside the commutator in Eq. (19c) is a matrix 
function only of A and r, and hence commutes with A. Substituting pr = +Pr into the term 
multiphed by r in Eq. (19b), and using TiCp^— p^TrC — 0, the traceless part of Eq. (19b) 
reduces to 

= {-rp'jrCp'^ + \pr,qs]+C6rs)AV ■ (20a) 
Since {C)av — ^eff^) this equation can be rewritten as the exact relation 

{[Qs:Pr])AV = ieS^^rs - T (^'jltCp'^) AV , (206) 

showing that the ensemble averaged commutator of the canonical coordinate and momentum 
operators has the form of the usual quantum mechanical canonical commutator, with i^s 
playing the role of the imaginary unit, and with a correction term proportional to r that is 
quadratic in the non-classical parts of the canonical momenta. Using Eqs. (14b) and (16c), 
Eq. (20b) can also be written in the form 

^{[QuPuDaV {Straus - StJrs) = -r(p's'TtCp'^)AV ■ (20c) 
t,u 

The derivation given here sharpens that given in Adler and Millard [4] , both in that here no 
approximations have been made, and that because we are working in complex Hilbert space, 
we have not had to first project out the parts of Qs and Pr that commute with ieff-* 

We conclude this section with several remarks on the principal result of Eq. (20b). 
First of all, the fluctuations about the ensemble average are fundamental to the possibihty 
of emergent quantum behavior. The average of the commutator on the left hand side of 
* As discussed in the Appendix, this projection is needed in real and quaternionic 
Hilbert space. 
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Eq. (20b) is not the same as the commutator of the averages, which in fact vanishes since 
{(1s)av and {pr)Av are both functions only of r and A, 

{[<ls,Pr])AV 7^ [{(Is)aV, (Pr-)AV] = . (21a) 

Secondly, in order to have emergent quantum behavior, the dynamics must be such that the 
second term on the right hand side of Eq. (20b) is much smaller than the first term on the 
right hand side, that is, one must have 

h-^T\l^',TiCp',)Av\«l ■ (216) 

We shall now show that this condition cannot be satisfied if C is an extensive thermodynamic 
quantity that grows linearly with the size of the system. To see this, we note that a second 
Ward identity, similar in form to Eq. (20b), can be derived by starting from 







J dfxSp^ [e-^«-^^^Trieff^p.] , (22a) 



in which the factor of C multiplying pr has been replaced by its ensemble average iesh, and 
then proceeding as in Eqs. (17b-20b) above. The resulting Ward identity is 

0^iMs-r{p',TrUshp'^)Av , (226) 

and is an analog in our context of the usual equipartition theorem of classical statistical 
mechanics. Now if C were an extensive quantity, the difference between C and its ensemble 
average if^gh would be a fluctuation that vanishes as the system size becomes infinite, which 
would make it permissible to accurately approximate the right hand side of Eq. (20b) by the 
right hand side of Eq. (22b). This would lead to the conclusion {[qs,Pr])Av = 0, that is, the 
thermodynamics would give emergent classical, rather than quantum mechanical, behavior. 
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Thus, the inequahty of Eq. (21b) can hold only if C is not extensive; this conclusion 
is consistent with the observation that to have the average of C play the role of the intensive 
quantity iefr^, one would expect that C should behave as a thermodynamically intensive 
quantity. We shall give further evidence for this conclusion in the next section, where we 
consider systems containing fermions as well as bosons. 

Although we have first discussed the purely bosonic model of this section for expos- 
itory reasons, it is easy to see from Eqs. (20b, c) that the inequality of Eq. (21b) cannot 
hold in a purely bosonic system. When the right hand side of Eq. (20c) can be neglected, 
multiplying by 5rs and summing over r, s gives the relation 

Y,^[quPu])Av5tu{l-N)^Q , (22c) 

t,u 

which ior N > 1 implies that {[qs-iPr])Av — 0, again giving classical behavior. Thus, a purely 
bosonic matrix dynamics system in complex Hilbert space cannot have emergent quantum 
behavior. However, we shall see in the next section that in the interesting case in which 
the numbers of bosonic and fermionic degrees of freedom are equal, a condition that holds 
for supersymmetric theories, the relation of Eq. (22c) is modified, and emergent quantum 
behavior becomes possible. 

4. Corrections to the Full Canonical Algebra in a 
General Unitairy Invariant Matrix Model with Classical Sources 

Although one can give derivations similar to that of the previous section for other 
canonical commutators (e.g., the [pr,Ps] and [gr,?s] commutators, or the [ps, ?r-] commu- 
tator obtained by interchanging the roles of p and q in the above derivation), and their 
fermionic anticommutator analogs, it is most efficient in the general case to use the sym- 
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plectic formalism introduced by Adler and Millard in Ref. [4]. In this notation one defines 

xi = qi,X2 = Pi,X3 = q2,X4 = p2, ...,X2R-i = Qr, X2R = Pr for the matrix phase space 
variables; in terms of these, the Hamilton equations of Eq. (5b), the generahzed Poisson 
bracket of Eq. (6a), and the conserved operator C of Eq. (7) take the form 

. _^ m 

Xr — UJrs r. 
r,s=\ 

— ^^^^ OCfLUfgXg , 
r,s 

with uj a numerical symplectic metric given (in terms of standard Pauli matrices Ti^2,z) by 
iT2 for a bosonic pair of canonical variables, and by — ti for a fermionic canonical pair (see 
Eqs. (lOa-c) on p. 202, and Eqs. (lOa'-c') on p. 224, of Ref. [4].) This symplectic metric 
obeys [4] the useful identities 

(236) 



^rs^rt — iOgrCOtr — ^st ■ 
r r 

We shall now consider a matrix dynamics generated by a general trace Hamiltonian 
H, that can contain fermionic as well as bosonic degrees of freedom, and in the statistical 
partition function shall allow the presence of nonvanishing classical sources for the classical 
parts of the phase space variables [cf. Eqs. (12a-c) above.] Thus we start from the ensemble 

p ^Z-^ exp(-TH - TrAC - J] JX) 



Z=J rf/xexp(-rH-TrAC'- ^ J, 



(24) 



which is now used in place of the ensemble of Eq. (13) in the definition of {0)av given in 
Eq. (15). For this ensemble, we consider the Ward identity obtained by using shift invariance 
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of the integration measure d/i starting from 

Z{TrCatx[)Av^ J diie'^^-^'^^-^r^r^rT^Catx'^, (25a) 

where at are a set of c-number auxihary parameters that are complex for bosonic Xt, and 
complex Grassmann for fermionic Xt- As in Eqs. (12a, b) above, we use the notation x[ to 
denote the noncommutative part of Xt that remains when the classical part is subtracted 
away, that is, x'^. = Xt — x^. 

The Ward identity derivation now proceeds exactly as in Eqs. (17b-20b). Making a 
constant shift of the noncommutative part x'g of the phase space variable Xg-, we have 
0=y" [e--^-^~^^-T.rJr<r^Catx'^ 
^ j d/,e-^H-TVAC-E.^.=x? (256) 

X [(-T^^'^H - TA5^,C)i:YCatx't + TY{5^,C)atx[ + TxCatSstSx'^ 
Now from Eqs. (23a, b) we have 

H =Tr (^^^ 5x'^ ^ ^ x'^^rsK , 

Tl\5x'fi =1l\^'^UJrs{x'r5Xg — erSx^X^) = Tr[A, ^ CUr-jX^j^X^ , (26) 
r r 

Tr{Sj:'^C)atXt^Tr[atXt,^uJrsx'^]Sx'g . 

r 

Substituting these into Eq. (25b) and multiplying by Z~^, we get 

={{—TTr'^2 K'^rs^x's ~ Tr[A, '^2'^rsx'r]Sx'g)TrCatXt 

' ~ ' (27a) 

+Tr[atx[, 22 ^rsx'rl^x'g + TrCatSst5x'g)Av , 

r 

which since Sx'g is an arbitrary traceless matrix (with the adjointness properties of x'^), 
implies that the traceless part of the operator multiplying 6x'g inside the trace must vanish, 

XI ^'r^''^ ~ XI ^rsXr])TrCatx[ 

+ [atx[, 22 ^rsx'r] + Cat6st)Av ■ 

r 
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Since A is a constant matrix, as before it can be taken outside the ensemble average, and so 
the second term in Eq. (27b) takes the form 

-[\,Y,^rs{x'rT^Catx',)Av] , (27c) 

r 

which again vanishes since the ensemble average inside the commutator is a matrix function 
only of A and r. Contracting the remainder of Eq. (27b) with '^^cuus, using Eqs. (16c) and 
(23b), and noting that [x'^, atx[] — [x^, crfXt] because the classical parts do not contribute to 
the commutator, we get as our final result 

{[xu, crtXt\)Av = ieS^utCTt - T {x'^Tr Catx[)Av ■ (28a) 

Equation (28a), like Eq. (20b) of the preceding section that it generahzes, is exact. 

In order to have emergent quantum behavior, it is necessary that the second term 
on the right hand side of Eq. (28a) be much smaller than the first term, that is, we require 

h-\\{x'^TrCatx'^)Av\«'^ ■ (286) 

Again, by replacing C by its expectation value at the start of the derivation leading to 
Eq. (28a), we get a second Ward identity 

= i^shujutc^t - T{x'^Txi^shatx[)Av ■ (28c) 

Hence the inequality of Eq. (28b) can be satisfied only if C is not an extensive quantity, 
since if C were extensive one could, in the large system limit, approximate it in Eq. (28b) 
by its expectation iesfi-i giving an expression that, by Eq. (28c), cannot be small. 

Because the derivation of this section is valid for fermions as well as bosons, one 
can in fact make a stronger statement about the conditions for emergent quantum behavior. 
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Letting the indices t and u in Eq. (28a) be either both bosonic or both fermionic, we get the 
respective relations 

{[Qr,Pr])Av ^iesfi " T {q[.TrCp'j.) Av r bosonic 

(29a) 

{{Qr:Pr})Av ^iesh - T {q'^T^C p'^) Av f fermonic . 
Substituting this into Eq. (7) for C, taking the ensemble average, and using Eq. (16c), we 
get 

ieflf^ =(C')Ay = ( Y [(lr,Pr\- ^ {^r,Pr})AV 
r bosons r fermions 

X / A ~ ^^^^^ 

E " E ) ^^ff^ ~ ( E ~ Yl ) {&Cp'^)av , 

\r bosons r fermions/ \r bosons r fermions/ 

which on division by h and transposition of terms gives 

E - E ) f^-'^{^rT^Cp'^)AV = ^eff ( E - E -1 ) ■ (29c) 

\r bosons r fermions/ \r bosons r fermions / 

When the condition of Eq. (28b) for emergent canonical behavior is satisfied, the left hand 
side of Eq. (29c) is a sum of very small terms. Assuming that this sum yields at most a finite, 
bounded total, let us consider the case in which r includes the spatial label of a translation 
invariant field theory. Then the number of bosonic and fermionic modes per unit volume 
contributing on the right hand side of Eq. (29c) must be equal, since if not, the right hand 
side of Eq. (29c) would become infinite as the spatial volume grows to infinity, contradicting 
the boundedness of the left hand side. Therefore, in a complex Hilbert space*, a candidate 
pre-quantum mechanics theory must have equal numbers of bosonic and fermionic degrees of 
freedom, making it plausible that such a candidate theory should be super symmetric. When 



* For a discussion of how our arguments must be modified in real and quaternionic 
Hilbert space, see the Appendix. 
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the numbers of bosonic and fermionic modes are in balance, Eq. (29c) simplifies to 

i J2 - Yl ]fi~'r{q',TrCp',)Av^-ieS , (29d) 

\r bosons r fermions/ 

showing that the remainder terms in Eq. (28a), that are neglected when Eq. (28b) is satisfied, 
sum in Eq. (29c) to give a total of unit magnitude. 

Corresponding to the Ward identity of Eq. (28a), we can derive a class of more 
general Ward identities by replacing atx[ in Eq. (25b) by a general U, constructed as a Weyl 
ordered (i.e., symmetrized) polynomial in the products {arXr}, with coefficients that are 
c-number functions of 1 and of ies- In place of Eq. (27b), we now get 

^ ~ (30a) 

r each Xs in U 

As long as U has coefficients that depend only on the matrices (or operators) 1 and ies, 
the second term on the right in Eq. (30a), which involves a commutator with A, vanishes 
by the same arguments as before. Contracting the remainder with ujus, but making no 
approximations, we get as the exact general Ward identity analogous to Eq. (28a), 

{[xu,U])av^{J2^us Y1 U{oneXs^Cy)Av-r{x',TrCU)Av ■ (306) 

s each Xs in U 

Suppose now that we can make the following two approximations, (i) we replace C in the 
first term on the right hand side of Eq. (30b) by its ensemble average iesh, and (ii) we neglect 
the T term in Eq. (30b). We then are left with the relation 

{[xu,U])av = (^(^us Y [/(one h'))Av , (30c) 

s each Xs in U 

which extends the effective canonical algebra inside ensemble averages to the commutator of 
Xu with a general Weyl ordered polynomial U. 
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By the methods of Appendix E of [4], Equation (30c) can be extended to include 
sources for the remainder parts x'^ of the phase space variables. Specializing the relations ob- 
tained this way to the case U — iefr implies that iefr can be freely commuted with phase space 
variables inside ensemble averages. As argued in [4], the resulting set of Ward identities then 
yields the canonical generator structure, including the time evolution relations, of Heisenberg 
picture quantum field theory. The only assertion in [4] that cannot be derived this way is 
the claim that the time evolution equation is more exact than the other generator relations; 
this claim used the assumption that C can be replaced by its ensemble average inside the 
T term, which we have seen is not correct. The remainder of the conclusions of [4] rest on 
the two approximations that we made above, which can be rephrased as the assumptions 
that, (i) in the terms of Eq. (30b) that involve the unvaried canonical ensemble p of Eq. (13) 
with a factor of C in the integrand, the fact that the ensemble is sharply peaked around the 
mean allows us to replace the integrand factor C by its ensemble average iesfi, and (ii) the 
canonical ensemble displays a certain rigidity, in that terms of the form J d/iSpTrCU can be 
dropped. On the other hand, we have seen that terms of the form J d/iSpTrieshU cannot be 
dropped; this does not contradict our assumptions because 5p can be rapidly varying around 
the peak of the ensemble. 

5. Discussion 

We have seen that a statistical mechanics can be formulated for a wide class of 
matrix models with a global unitary invariance, and that within this statistical mechanics, 
the ensemble averages of canonical variables obey the exact relations of Eq. (20b), (28a), 
and (29a), and (30b). When these can be approximated by dropping the r terms [and, 
in Eq. (30b), replacing the C insertions in U' by their ensemble averages], the result is 
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emergent quantum mechanical behavior for the statistical ensemble averages. The condition 
for validity of the approximation of neglecting the r terms is rather delicate: we have argued 
that it requires that C should be an intensive thermodynamic quantity, and that the numbers 
of bosonic and fermionic degrees of freedom should balance. 

We conclude with some brief remarks: 

[1] In the first two references of [3], it is shown that one can readily formulate trace 
dynamics models in which global unitary operator invariance is gauged to give a 
local unitary operator invariance. Since global unitary invariance is a special case of 
local unitary invariance, the considerations of this paper apply to these models. 

[2] In Refs. [5], it is shown that supersymmetric Yang Mills theory, and the related 
"matrix model for M theory" , fit naturally into the trace dynamics framework ana- 
lyzed in this paper. In these models, C vanishes up to a surface term contribution, 
a behavior consistent with its being an intensive thermodynamic quantity. 

[3] Although our final results of Eqs. (20b), (28a), and (29a) superficially resemble the 
string-inspired formula of Eq. (lb), there is an important difference. In Eq. (lb) the 
coordinate is a quantum operator, as is usual in nonrelativistic quantum mechanics, 
whereas in our results of Eqs. (20b), etc., the coordinate is merely a degree of freedom 
label r, as it always is in quantum field theory, and the coordinates and momenta 
are canonical field variables with label r. It may be possible to make a connection 
between the two types of modified commutation relations when the metric structure 
of the coordinate manifold is taken into account, using the fact that the proper 
distance is related to the coordinate interval by ds^ = g^ydx^dx^ . In a field theoretic 
interpretation, dx^ is just the change in the degree of freedom label, but the metric is 



a dynamical variable, and hence so is ds. This suggests that there may be an analog 
of Eq. (lb) involving the relativistic proper distance, and that this is the relation to 
be compared with our results in this paper. 

[4] As is well known, in a complex Hilbert space the canonical algebra [q,p\ — i, [q, i] — 
\p,i] =0 cannot have finite dimensional representations, since this algebra implies, 
for example, the relation q^p^ -\-p^(f — 2qp'^q — —2, which in a finite dimensional 
Hilbert space would have a left hand side with trace zero and a right hand side with 
trace nonzero. However, it is consistent for the canonical algebra to emerge as the 
limit A?^ ^ oo of an algebra in an dimensional Hilbert space, which is the behavior 
argued for in Ref. [4] and here. Because the emergent canonical algebra involves 
not the imaginary unit i of the underlying complex Hilbert space, but rather the 
operator igff with Triefr = 0, a basis for the operator algebra in the emergent theory 
is provided by a set of operators that commute with ies, together with one additional 
operator that anticommutes with i^s, and that plays the role of the time reversal 
operator in the emergent complex quantum mechanics. In fact, because the condition 
Trieff = implies that one can find a representation in which i^s is a real matrix 
(just as for Pauli matrices pi,2,3 the matrix ips can be given the real form ip2 by a 
change of representation), the quantum mechanics emergent from matrix dynamics 
has the structure of a complexified real quantum mechanics, for which the operator 
algebra has the form just described (see, e.g.. Sec. 2.6 of the second citation in Ref. 



[5] We have seen that the emergence of quantum mechanical behavior from matrix 
model statistical mechanics requires a certain "rigidity" of the statistical ensemble. 
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It is easy to see [8, 9] that this rigidity is a sufficient condition for the canonical 
and microcanonical ensembles to give the same Ward identities, and hence the same 
emergent quantum behavior. The need for a rigid statistical ensemble in our context 
suggests a possible analogy with the concept of London rigidity in the theory of 

— * 

superconductivity [10]. In the presence of an applied vector potential A, the induced 
current density in a metal is given by 

In a normal metal the two terms on the right hand side of Eq. (31a) nearly cancel, 
leaving a small residual diamagnetism. However, in a superconductor the rigidity of 
the wave function leads to the vanishing of (p ), giving perfect diamagnetism and 
the Meissner effect. An analogy with the results of this paper would equate normal 
metal behavior with the case in which C can be replaced by its ensemble average in 
the T term; in this case the right hand side of Eq. (20b) is approximately equal to 
the right hand side of Eq. (22b), leading to cancellation of the emergent canonical 
commutator. Similarly, the analogy would equate superconducting behavior with the 
case in which the r term containing C can be dropped because of "rigidity" of Sp, 
leading through Eqs. (20b), (28a), and (29a) to an emergent canonical commutator 
as an analog of the superconductive Meissner effect. 
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Appendix: Real and Quaternionic Hilbert Space 

For reasons that we now describe, in real and quaternionic Hilbert space the argu- 
ments of this paper must be modified and yield weaker conclusions. The underlying reason 
for this modification is that only in complex Hilbert space can one have a non-real trace 
that nonetheless obeys the cyclic property. In quaternionic Hilbert space, as a consequence 
of the noncommutativity of the quaternions, only the real part of the ordinary trace obeys 
the cyclic property. In real Hilbert space, the trace is necessarily real, and the trace of any 
anti-self-adjoint operator vanishes. Thus in these two cases, if one follows [3, 4] and defines 
the graded trace Tr to be the one that obeys the cyclic property, then the definition must 
include taking the real part, and to get a nonzero result one must require the operator argu- 
ment V of TrV to be self-adjoint. As a consequence, in the general case derivation analogous 
to that starting from Eq. (25b), one must consider Trjieff , C}atX^ rather than TiCatX^. This 
gives the following analog of Eq. (28a), 

(K, CTtiieS: Xt}])AV = {^eflf, ieS^YuJutCTt " T{x'^Tr{i^s, C}atx[)AV , (^1) 

where the prime on the first term on the right hand side of Eq. (Al) indicates extraction of 
the traceless part. However, since the traceless part of {ies,iesfi} = — 2/i is zero, Eq. (Al) 
becomes 

{[xu:Crt{ieS:Xt}])Av ^ -r{x'yTr{ieS,C}atx[)Av , (^2) 

which has a structure more hke the equation obtained by subtracting Eq. (28c) from Eq. (28a) 
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than like Eq. (28a) itself. In complex Hilbert space, Eqs. (A2) and Eqs. (28a, c) all hold, 
and conditions for emergent quantum behavior can be formulated from the latter as in the 
text; in the real and quaternionic Hilbert space cases, only Eq. (A2) holds, and we cannot 
proceed with the analysis of the text. 

The remarks made here correct a subtle error made in Sec. 6 of Adler and Millard [4] . 
There, the classical sources introduced in this paper were set to zero, and Ward identities 
were derived by varying Xs, not just the noncommutative part x'g as we have done here. This 
gives as the Ward identity analogous to Eq. (Al) 

{[Xu, CTtiieS, Xt}])AV = {^eflf, ieS^^utC^t " T(i;„Tr{ieff, C}atXt) AV ■ (^3) 

Equation (A3) is a valid relation, but it actually implies two relations, quite different in 
structure, for its classical or c-number part and its primed or traceless part. Separating 
Eq. (A3) into these parts, and ignoring interference terms between x'^ and x' in the r term, 
one finds that the primed part of Eq. (A3) gives Eq. (A2), while the c-number part of 
Eq. (A3) becomes an equipartition identity for x'^ and gives no direct information about the 
expectation of the commutator. Hence the arguments of [4] do not, in their present form, 
provide evidence for emergent quantum behavior in the real and quaternionic Hilbert space 
cases. 
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